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li. Aiitfoct This paper describes improved methods of calculating 
wall interference corrections for the ONERA large windtunnels. 
The mathematical description of the model and its sting suppor.-ti 
have become more sophisticated. An increasing number of singurj 
larities is used until an agreement between theoretical and ex-j 
perimental signatures of the model and sting on the walls of the 
closed test section is obtained. ‘ The singularity decentering 
effects are calculated when the model reaches large angles of 
attack. The porosity factor cartography on the perforated wal]s 
deduced from the measured signatures now replaces the reference 
tests previously carried out in larger tunnels. The porosity 
factors obtained from the blockage terms (signatures at zero 
lift) and from the lift terms are in good agreement. In each 
case (model + sting + test section) , wall corrections are now 1 
determined, before the tests, as a function of the fundamental j 
parameters M, CS, CZ . During the wind tunnel tests, the cor- 
rections are quickly computed from these functions. i 
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ABSTRACT 


The methods used to compute wall interference corrections 
for the ONERA large windtunnels have been improved over the years. 
The mathematical description of the model and its sting support 
is more and more sophisticated; an increasing number of singular- 
ities is used until an agreement between theoretical and experi- 
mental signatures of the model and sting on the walls of the closed 
test section is obtained. The effect of the singularity displace- 
ment from the central position is calculated when the model reaches 
large angles of attack. 

The porosity factor cartography on the perforated walls de- 
deduced from the measured signatures now replaces the reference 
tests previously carried out in larger tunnels. The porosity 
factors obtained from the blockage terms (signatures at zero lift) 
and from the lift terms are in good agreement. 

In each case (model + sting + test section) wall corrections 
are now determined, before the tests, as a function of the funda- 
mental parameters M, CS, CZ . During the windtunnel tests, the 
corrections are quickly computed from these functions. 
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INTRODUCTION 
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The final closing discussion of the AGARD/FMP Symposium at 
Valloire in 1975 [1] confirmed the need to apply wall interfer- 
ence corrections to tests carried out in existing transonic wind- 
tunnels. The aim of these corrections is not to give the same 
results as those found during flight, because of the different 
Reynolds numbers, but to restore the speeds which have been mod- 
ified by the presence of walls, both in magnitude and in direction. 
In this way, the results obtained using the same Reynolds for var- 
ious windtunnels and for different model scales should be homogen- 
eous. These types of comparison are a precious guide for checking 
the calculations of wall effects. Numerous comparative programs 
have thus become more and more sophisticated: this was the case 
with the tests carried out on standard ONERA models [2] , or tests 
now being carried out on the F4 model within the framework of the 
GARTEUR group, made up of various organizations. Such is the con- 
stant concern of ONERA as it multiplies the comparisons of tests 
carried out in the various large three-dimensional Si, S2, S3MA 
and Fl windtunnels [3]. 

The investigation of wall effects has given rise to abundant 
literature since 1919. A study [4] of these effects on perforated 
walls has resulted in a list, although certainly incomplete, of 77 
published references in 1977. At least as many works on slotted 
walls have been published. This large number of publications would 
lead one to think that the problem is solved. It appears, however, 
judging from the list of congresses held on the subject over the 
past ten years, that this is not the case. The AGARD/FDP congres- 
ses held in London in 1975 [5], Rhode Ste Genese in 1976 [6], and 


*Numbers in the margin indicate pagination in the original test. 
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Munich in 1980 [7] presented 27 reports on wall effects [8]. The 
last AIAA conference held this year at Williamsburg had a full ses- 
sion on wall effects [8] . We could, of course, mention this con- 
gress where 17 reports on wall effects only, will be given. 

1 - RECENT IMPROVEMENTS IN WALL INTERFERENCE CALCULATIONS AT ONERA 

The work on adaptive wall studies currently being conducted 
in ONERA two-dimensional windtunnels will not be considered here. 
This report will be concerned only with three-dimensional models 
in standard cylindrical test sections used for industrial tests. 

The reasons for conducting investigations to improve wall inter- 
ference calculations in such cases may be classed in five groups: 

-modern aircraft projects require greater accuracy, particularly 
for civil transonic transport aircraft 
-larger angle of attack domains for modern fighters 
-more powerful computers available for calculations 
-measurements of signatures of model-sting combinations on the 
walls of test sections 

-increased model sizes versus test sections offering a better 
representation of certain details, and a higher Reynolds number 
for testing. 

Recent improvements in wall interference computations per- 
formed at ONERA will be described in this publication according 
to three viewpoints. 

A) Improvement of the mathematical descriptions of models , in 
terms of both the number of singularities, and their location, par- 
ticularly for cases of high angles of attack; actually because 
mechanisms are inserted in the test section to allov; for models 
v/ith higher angles of attack, these models become displaced from 
their central position. This report also describes the model 
supports, v/hereas, until now, models have been studied separately. 


B) Utilisation of measurements of signatures of model-sting 
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combinations on the walls of test sections in order to, first, 
check out the validity of the mathematical description of model- 
stinq combinations in the closed test section, and, second, to 
deduce the porosity factor cartographies on ventilated walls. This 
new method of procedure, therefore, makes it possible to omit the 
reference tests, such as testing in closed test sections, or in 
larger windtunnels, which have been used until now [9], 

C) Reviewing the methods used to compute wall interference 
corrections: the tables of correction coefficients calculated in 
advance for a given test section, as a function of the numerous 
parameters, such as wall porosity, relative model span, wing sweep, 
span lift distribution, etc. - will be discarded. The corrections 
are computed for each case (test section-model-sting) for a Mach 
grid, Cz, Cx, and the corrections are deduced in the form of poly- 
nomials as a function of these 3 parameters. This procedure offers 
a rapid method of computing the corrections for each test point. 

2 - MATHEMATICAL DESCRIPTION OF THE MODELS AND STINGS /3 

2.1 - Representing the Volumes 

The volimie of the model was initially [4] represented using 
a single doublet situated at the center of the model which is lo- 
cated iii the standard position at the quadrant of the mean aero- 
dynamic chord. Empirical coefficients were added to account for 
the fuselage slenderness, and the relative dimensions of the model 
in the test section. This made dt possible to calculate the Mach 
number, normally assumed to be at the center of the model. Con- 
versely, the volume blockage gradient is considered to be incor- 
rect, where one singularity cannot account for the law of area. 
Moreover, the model signatures on the walls of the test section 
were too concentrated, and this led to the estimation of a pre- 
mature blockage in the closed test section. After equipping the 
walls of the test sections v;ith pressure pick-ups, it was found 
that the experimental signatures were noticeably different than 
the theoretical signatures. 
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For these various reasons, the fuselage of the models was 
represented by a set of doublets, the fuselage length is divided 
into N equal intervals for which the elementary volumes are cal- 
culated from the law of area together with their respective center 
of gravity. Each elementary volume is represented by a doublet 
situated at the center of gravity of the volume; the doublet in- 
tensity, proportional to the elementary volume, is assigned a term 
of compressiblity deduced from the corrected Mach number calculated 
in its position after several iterations. The doublets number N 
retained is that above which the signatures and corrected Mach dis- 
tributions on the fuselage axis are not modified. 

Figure 1 shows the signatures of a slenderness ellipsoid, _in 
a closed test section, obtained with a description by 1, 2, 3, and 
12 doublets. The signature becomes asymptotic after 12 doublets. 
After this number, and especially with a single doublet, the sig- 
natures are too concentrated and do not exhibit the stat of spread 
obtained v;ith a large enough number of doublets. 

A simple rule, deduced from calculations for various ellip- 
soid slendernesses, consists of taking a number of equal or double 
doublets of the fuselage slenderness. 

Figure 2 gives a comparison of the calculated and measured 
signatures of a large scale missile 5.3 meters long, obtained in 
a closed test section of the SIMA windtunnel 8 meters in diameter. 

In this case, the satisfactory agreement between these signatures 
is obtained v/ith 30 doublets. The calculated signatures "custom- 
ize" the model configuration quite well, as shown by the measure- 
ments given in the case of a single fuselage or of the full missile 
with nacelles and tail units. The signatures reflect the form of 
the law of area. The signatures calculated with a single doublet, 
given for reference, in a broken line, would be perfectly incor- 
rect, because they are too concentrated and have a, maximum Mach that 
is too high. 




Fig. 1 - Effect of the number 
of doublets on the signature 
of an ellipsoid in a closed 
test section. 


theor. 30 doubMts 
» 1 • 



Pig. 2 - Theoretical and experi 
mental signatures of a missile 
in a closed test section. 


From this fact, the blockage limits of the model in the test 
section were obtained prematurely by the calculation (figure 3), 
whereas with a large enough number of doublets, the blockage is 
predicted at a Mach number exceeding about 0.02. The tests sub- 
stantiate fairly well that the blockage appears only above Mach 
0.9 for the complete missile. 


The distribution of the corrected Mach number along the model 
axis is directly proportional to the local coefficients 523 ( 4 ] for 
a single doublet or v/ith analogous coef f icents* asz for N doublets. 
Figure 4 gives, for different factors of reduced porosity ^ 

the longitudinal distributions of coefficients iii and 'iilsc 'for 1 

> • i 



and 20 doublets representing the 
slenderness ellipsoid 6 with 1% 
obstruction. With 20 doublets, 
the curves are always smaller. 

In this same figure, the Mach 
corrections (Me - Mu) , at the 
ellipsoid center, are also re- 
duced by using the correct num- 
ber of doublets. 

It follows that the "Arch- 
imedean thrust" corrections, due 
to volume blockages , themselves 
directly related to the longitu- 
dinal gradients of the coef- 
ficients, i.e. 'a's' andl^'sc for 1 
and N doublets, are modified (fig. 
5) . These corrections are 

therefore reduced to an abso- 
lute value (figure 5) . 

If we intuitively situated the previously considered single 
doublet at the center of this volume for the case of an ellipsoid, 
there would still be an unkown when a single doublet is used for 
a fuselage; its location. From this standpoint, the correc- 

tion changed (figure 6) , v;hereas with a set of the correct number 
of doublets, this correction stayed the same. 

2.2 - Representing the Wake and the Separations 

For the case of a fuselage, the v;ake blockage was calculated 
by using a source, situated at the center of the model, of intensity 
proportional to the theoretical drag coefficient. In recent calcu- 
lations, the source may be placed at any location. The effect of 
the backward travel of the source from the center of the model at 
the base is shov/n in figure 7: it reduces the Mach corrections 
along the fuselage axis together with the level of the signatures 


tftaor. 1 doublet \ 

___ tbeor. 30 , 

• ^ exper. 



Fig. 3 - Prediction of the Mach 
blockage number in the closed 
test section. 
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Fig. 4 - Influence of the modeling of an 
ellipsoid on the Mach correction on ventil- 
lated walls. 



Fig. 5 - Influence of the modeling of an ellip 
soid on the buoyancy correction on ventilated 
walls. 
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Fig. 6 - Influence of the doublet 
position on the buoyancy correction. 
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on the wall of the test 
section. It is there- 
fore important to specify 
the location of the wake. 

For the case of an 
aircraft model, the pro- 
blem of the location and 
of the singularities used 
to represent the wake and 
the separations with large 
angles of attack is defin- 
itely more complex. The 
best models today have 
been proposed by HACKETT [9, 

10 ] . 
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Fig. 7 - Influence of the posi- 
tion of the source and the cal- 
culation iterations on the sig- 
natures. 


It should be pointed out, as 
already mentioned in section 2.1, 
that the singularities used to 
represent the volumes and the 
wakes have weighted intensities by 
the compressibility term deduced 
from the Mach number corrected at 
the location of each singularity. 
This implies an iterative calcula- 
tion: the number of iteration, of 
course, increases with the Mach 
number. Figure 7 illustrates the 
effect of the iterations for the 
case of the preceding missile at 
SIMA at Mach 0.90. The importance 
of these iterations is obvious on 
the calculation of blockage in the 
closed test section. 



2.3 - Inclusion of the Sting Support 
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Until now, the model has been described mathematically with- 
out its sting support. Figure 8 shows that for the case of a 
standard ONERA MS model [2] mounted on a straight sting in the 

transonic test section with 
closed configuration of the 
S3MA at Mach 0.83, the theo- 
retical signatures at zero 
lift are different than the 
experimental signatures using 
20 doublets to represent only 
the model. The difference 
becomes more significant to- 
ward downstream. 

The modeling of the sup- 
port using 15 doublets added 
to that of the model, made it 
possible to obtain an excel- 
lent agreement between the 
theoretical and experimental 
signatures, even outside the 
area occupied by the model. 

It is therefore important to 
account for the area requirement of the sting support, particularly 
because the sting areas often overshoot the model's partition. 

The corrected Mach distributions, along the axis of the model, 
and proportional to the Siax , coefficients, are given in figure 9. 
These are comparative effects of the interaction terms due to the 
wall effects (and not effects of the set of potentials due to the 
model being in an unlimited atmoshpere and due to the ' ‘I’c ' inter- 
actions. The blockage effects from the sting, downstream from 
the model, v/ill have repercussions here in the entire area occupied 
by the model. This influence increases as the Mach decreases. 


S3Ma W.X(clos3d walls) M.0.S3 



Fig. 8 - Model and Sting Signatures 
in a closed Test Section. 
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Fig. 9 - sting Effects on the Blockage Corrections 

The derivative ^'sz > of the preceding function s-sxtx)'. leads 
to the correction of the Archimedean thrust._AC«_ induced by the 
blockage terms. Figure 9 shows the alternation of the y j . 

curves. For the model only, this curve had an odd number, which 
leads to a correction "acxs of practically zero. Conversely, the 
presence of the sting gives which are almost always positive, 

except on the tail region. Consequently, the ac<s correction is 
negative. The resulting error, by not taking the sting into con- 

-4 

sideration, is about 10.10 in drag, and is therfore large. The 
influence of the sting support on the model in an unlimited atmos- 
phere potential) should be added to this correction. 

2.4 - Representing the Lift 

The lift is represented by a vortex sheet, without inclina- 
tion, by accounting for the relative span, its sweep and the lift 
distribution on the span, if it is knov/n [4]. Such a representa- 
tion is appropriate for moderate angb'e.s of attack, but a few pre- 
cautions should be taken for large angles of attack [3] . 

This modeling may be checked by noting the differences in the 
signatures on the high and low wall of the test section. Indeed, 
a comparison of the slopes on a point X, Y of the test section of 



the curves showing the differences beween local Machs measured on 
the 2 slopes M„ M_, as a function of C , is obtained using a S'* , 

no Z 

coefficient analogous to the angle of attack correction coefficient 
'-T; [4]. 

Using the same standard ONERA M2 model/ figure 10 gives a 
comparison of the theoretical and experimental tSr/ coefficients 
which coincide satisfactorily for the cartography of the signatures 
on the two walls (curves at Y = 0 in the vertical, plane of symmetry 
and at X = 0 in the transversal plane) . This validates the model- 
ing of the lift terms. 

3 - DETERMINATION OF THE POROSITY FACTORS ON THE WALLS 
3.1 - Signatures At Zero Lift 

The mathematical representations of the models and stings 
are now sophisticated enough to obtain valid signatures in the 
closed test section. The problem of the boundary conditions on 
perforated walls, therefore, may now be handled. 

Currently, the direct method which would consist of directly 
deducing the wall interference corrections from the theoretcial 
signatures is used only for 3D flow. An indirect method, by means 
of the porosity factor cartography on the walls, is the preferred 
method at the present time, because a basic objective searched for 
is to show whether the concept of uniform porosity should be aban- 
doned or not. 

The signatures are calculated as a function of a constant /8 

porosity factor on the walls. The theoretical signatures are 
reported on such a signature grid as a function of X and Q. Fig- 
ure 11 thus gives the experimental signature, at zero lift, for 
the same ONERA M2 standard model mounted on a straight sting in 
a closed test section and in a test section with horizontal per- 
forated walls. It may be seen that the latter signature coincides 
remarkably v;ell with the experimental signature for Q = 0.2 which. 
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Fig. 10 - Comparison of Lift Terms Effects 
on the Signatures in a Closed Test Section. 


in this case, validates the concept of uniform porosity. Since /8 
this model is centered, the signatures used here at zero lift do 


not depend on the blockage terms. 
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Fig. 11 - Determination of the 
Porosity Factors From the Sig- 
natures on Perforated Walls 
(Blockage Terms) . 


3.2 Signatures At Nonzero Lift 

As indicated in paragraph 
2.4, the difference between 
the signatures on the 2 high 
and low walls may be used to 
check the validity of the 
lift modelings. Once the 
validity is confirmed in the 
closed test section, it may 
be applied to porous walls 
to obtain the porosity fac- 
tor cartography of the walls. 
Figure 12 illustrates this 
method for the case of the 
standard ONERA M2 model at 
S3MA in the perforated v/all 
configuration. 
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Fig. 12 - Determination of the Porosity Factors 
From the Signatures on Perforated Walls (Lift Terms) . 

In terms of the grid calculated as a function of the 

porosity parameter 9 it is deduced that in the vertical plane of 
symmetry/ a constant value of 0.2 is obtained for a which agrees 
with the value found at zero lift (figure 11). It may therefore 
be noted that there is an identity of the porosity factors of the 
blockage (Cz zero) and lift terms. On the span, there is a slight 
increase of Q from 0.2 to 0.35. 

It now appears that the porosity factor cartographies may 
be deduced from signatures on walls and that the determination of 
porosity factors in reference tests performed in a closed test 
section or in a larger windtunnel may be abandoned. This new 
method has tv;o incontestable advantages: first, supplementary 
reference tests are not needed, and second, the porosity re- 
quirements for every Mach and Reynolds range are defined, and 
this is often not possible during reference testing. However, 
this method requires a reading of the signatures which extends 
the testing time, although they can be limited to a number of Mach 
and angle of incidence points. 




,1 



Figure 13 - Quick Method for 
Calculating the 
Mach Number. 
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3.3 - Rapid Method For Determining the Corrected Mach 
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Based on the formulations of the theoretical signatures on 
the model on the walls [4] , an interesting method for determining 
the corrected Mach number has been obtained. This method consists 
of deducing the corrected Mach Me from the theoretical Mach Mm on 
the entire wall perpendicular to the strut of the model (along the 
wing extension) , whereas normally Me is calculated from the theo- 
retical Mach as far as possible upstream from the model to be free 
of its influence. In this case, however. Me is calculated from 
the most interacting point of the model. 

The incontestable advantage, as shown in figure 13, is to 
be able to do without knowledge of the wall porosity, because the 
correspondence between Me and Mm is practically independent from 
the porosity factor on this particular point. 

Only knowledge of the model's volume is needed. This rapid 
method requires only a very short computer time and makes it pos- 
sible to conduct a test at a constant Me. However, it requires a 
good modeling of the volumes of the models by taking into consider- 
ation the area requirements. Figure 13 thus shows that some 20 
doublets are needed for forming a model of the mock-up. This fig- 
ure illustrates the very small variation of the Me - Mm correction 
as a function of the porosity factor for an ONERA M5 model in the 
S2MA windtunnel. 

4 - THE EFFECT OF THE MODEL DISPLACEMENT FROM THE CENTRAL POSITION 

In the mathematical representations described in paragraph 2, 
all of the singularities used are situated along the horizontal 
median plane of the test section. Such modelings are applicable 
to the case of models which are centered in the test section and 
for moderate angles of attack. In this particular case, the block- 
age terms result from the volume and wake singularities, and the 
angular flov/ corrections are deduced from the lift singularities. 



If we now take models offset from their central position into 
consideration, or even models centered at a large angle of attack, 
we would have to take into consideration a set of singularities 
for any position in the test section, particularly outside the /lO 

horizontal median plane of the test section. The calculation of 
of the wall effects becomes complicated by quadratic interactions: 
accordingly, angular corrections are induced by volume and wake 
singularities and blockage corrections are deduced from the lift 
singularities. The importance of these terms increases with the 
singularity displacement from their central position; their effects 
are not negligible, and this is shown in the two examples given here. 



Fig. 14 - Influence of the Displacement From 
the Central Position of a Model on the Cor- 
rections in the Closed Test Section. 

The first example refers to the new ONERA fl windtunnel 
having a closed test section in the low speed range (figure 14) . 

Two! standard models, civil and fighter aircraft, of standard Fl 
size are examined. The Mach corrections are given as a function 

of the relative displacement from the center t/K of the models by 
taking into account, first, the volime and wake singularities, and, 
second, by adding the lift singularities; the significance of the 
latter terms is obvious. Figure 14 shov;s the lift gradient modi- 
fications resulting from the model's displacement from a central 




position. The lift curve spread before corrections depends on 
the type of model: the order of the curves on an increasing gra- 
dient would therefore be: negative offset, zero is positive for 
civil aircraft, zero offset, negative then positive for the fight- 
er. 


The second example refers to a civil aircraft model in the 
S2MA transonic windtunnel with horizontal perforated walls & a re- 
duced porosity factor q close to 0.7. Figure 15 gives the Mach 



Fig. 15 -Influence of the Offsetting of a Model 
on the Corrections on Perforated Walls. 


as a function of the offsetting at various test Mach, by including 
the wake and blockage singularities only, then by including the 
lift singularities; the curve dissymmetry increases with the Mach. 
The same applies to alterations of the lift gradient corrections. 


A test result obtained on the standard ONERA M2 model in the 
S3MA windtunnel with perforated walls (figure 16) , for 2 positions 
of the model; .centered and positive decentering, shows that besides 
the variations of the lift gradients and the stabilities, the de- 
centering modifies the angle of attack and the pitching moment at 



zero lift. For the centered models, the classical calculations 
of corrections cannot explain such deviations. Conversely, by 
calculating with the decentered singularities makes it possible 
to predict; an angle of attack correction at zero Cz is deduced 
from decentered volume and wake singularities: figure 17 shows 

I 

that this type of correction rapidly increases with the decenter- 
ing and the Mach. The computer program for corrections with de- 

I 

* centered singularities, in large enough numbers, is still being 

prepared. 

5 - CALCULATIONS OF WALL INTERFERENCE CORRECTIONS /II 

Owing to the increasing complexity of the calculations of 
wall interference corrections, resulting from more sophisticated 
modelings of the model-sting combinations, particularly for cases 
of decentering, it is no longer possible to construct tables of 
correction coefficients as done previously. 

^ The corrections are computed for each test section, model, 

V- sting support, before a test: they require more and more soph- 

e> isticated computer programmes. The corrections are defined for a 
grid of M, CX, CZ covering the test program and they are restored 
in the form of polynomials as a function of these three parameters. 
While the tests are being carried out, the corrections are computed 
for each test point from the previous polynomial, and are therefore 
clearly reduced and applicable to industrial testing programs. 

CONCLUSIONS; FUTURE PROSPECTS 

Figure 18 summarizes current improvements and those projected 
for the future. Regarding the mathematical modelings of the model- 
sting combinations placed in a test section, the transversal area 
requirements must still be defined to avoid concentrations of wing 
volumes, for example, on the model axis- Modelings of apex vorti- 
cities, nacelle jets and struts remain to be defined, as well as 
the means for deflecting the vortex sheets. It should be pointed 






Fig. 16 - Test Results on Models With 2 Decenter- 
ings on Perforated Walls. 





2h 


S2Ms «.Z 

Q.a? 

2b 


out that only the far fields are 
taken into consideration, and that 
consequently, only rudimentary, 
yet adequate, modelings are to be 
considered. 
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The indirect method of comput- 
ing corrections, with determination 
of porosity factor cartography, will 
be discarded and replaced by a dir- 
ect computer method using signatures, 
such as the one already used indus- 
trially for two-dimensional flow [11] 
This method is based on the measure- 
ment of pressures on the walls of 
test sections, and the minimum num- 
ber of pressures will first be de- 
termined for each type of correction in order to minimize the use 
of equipment on the walls, and therefore the data volume to be 
processed, and as a result, the length and cost of testing. 


Fig. 17 - Inlfuence of the 
Decentering of a Model on 
the Angle of Attack at 
Zero Lift. 
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Figure 18 - Conclusions: Results Acquired, Future Studies. 

The constant concern of windtunnel managers will result in 
multiplying the checks made of the corrections, by having recourse 
to every possible comparison of the results, v;ithin a given test 
section for various possible configurations, or between various 


windtunnels. 
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